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Sources of photon pairs based on the spontaneous parametric down conversion process are commonly used for 
long distance quantum communication. The key feature for improving the range of transmission is engineering 
their spectral properties. Following two experimental papers [Opt. Lett., 38, 697 (2013)] and [Opt. Lett., 39, 
1481 (2014)] we analytically and numerically analyze the characteristics of a source. It is based on a /3 barium 
borate (BBO) crystal cut for type II phase matching at the degenerated frequencies 755 nm —)► 1550 nm + 
1550 nm. Our analysis shows a way for full control of spectral correlation within a fiber-coupled photon pair 
simultaneously with optimal brightness. @ 2015 Optical Society of America 
OCIS codes: 190.4410,300.6190,270.4180,270.5565 


Sources of polarization entangled photon pairs based 
on the process of spontaneous parametric down conver¬ 
sion (SPDC) [l]-[^ are now commonly used in various 
experiments such as those testing fundamentals quan¬ 
tum mechanics, for quantum communication protocols 
or quantum information processing. Each application de¬ 
fines a set of specific requirements for the characteristics 
of photons for the best performance. For example, for a 
heralded qubit encoded in a polarization state any cor¬ 
relation in any other degree of freedom, within a pho¬ 
ton pair, will reduce the purity of the resulting state. 
Therefore there is a need to engineer characteristics of a 
source to produce uncorrelated pairs [£-11 . In turn, the 
polarization entanglement is a basic block of quantum 


communication 12 13] and quantum metrology 14 15 


The photons are transmitted through optical systems 
and single-mode fibers (SMF). Optical elements are dis¬ 
persive and cause unwanted modification of the photons’ 
temporal characteristics. The deteriorating dispersion ef¬ 
fects can be reduced by using photon pairs that are spec¬ 
trally narrow and feature appropriate spectral correla¬ 
tions 


16 18 


The spectra of SMF coupled photon pairs generated 
in SPDC process can be negatively 19 20 , positively 


16 18,21-23 or not correlated at all 24-26 . The ef¬ 


fect of negative correlation originates from the energy 
conservation relation and is a consequence of pumping 
by a continuous wave (CW) laser. Uncorrelated spec¬ 
tra can be achieved in a nonlinear crystal by using a 
pulsed laser pump with careful setup parameter choice 
or by using narrow spectral filters. However, for positive 
correlation, an additional requirement of group velocity 


matching needs to be fulfilled 21 


Here we perform numerical analysis of the SPDC pho¬ 
ton pair source experimentally analyzed in Refs m\ 


which is schematically depicted in Fig. We discuss 
trade-offs between the correlation parameter, pair pro¬ 
duction rate and SMF coupling efficiency, which are fig¬ 
ures of merit allowing to estimate source’s suitability for 



Fig. 1. A single photon source is based on a ;d-barium 
borate (BBO) nonlinear crystal of thickness L cut for 
type H phase matching at frequencies 755 nm ^ 1550 
nm + 1550 nm and half cone opening angle 3 deg. A 
pump pulse of Tp duration is centered around 775 nm 
and its transverse spatial mode is assumed to be a Gaus¬ 
sian function with a characteristic diameter 2wp. The 
photon pairs are coupled into single mode fibers (SMF), 
which together with coupling stages (C1,C2) define a col¬ 
lected transverse spatial modes, which are assumed to be 
a Gaussian functions of diameters 2wf. 


a given application. In the setup a femtosecond laser 
pulses centered at 775 nm are used to pump a BBO crys¬ 
tal, which optical axis orientation is set at 29.68 degree 
allowing type H phase matching at degenerate frequen¬ 
cies 1550 nm. The two created photons propagating at 
an 3 degree angle with respect to pump beam propaga¬ 
tion direction are coupled into SMFs. 

Firstly, the goal is to derive a simplified compact ana¬ 
lytical formula for characteristics of our source. Follow¬ 
ing the mathematical model presented in Ref. [Tq|[^ 
we use an effective phase matching function (EPMF) 
&{uJs,uJi) in order to fully describe joint effect of the 
pump beam spatial profile A^^^^(fc^), the crystal’s phase 
matching smc{AkzL/2) and the collected modes’ pro- 
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files Us{ks±), Ui{ki±). It is defined as the following over- 
lap; 0(ws,Wj) = / dksj_dki±Us{ks±)ui{kij_)Ag^°'\ks± + 
ki±)smc{AkzL/2). Here Akz = kpz — kgz — kiz stands 
for a phase mismatch of pump and resulting pho¬ 
tons. Then a state of a photon pair propagating in 
single mode fibers can be conveniently expressed as 
a product of a pump pulse spectral mode Ap{uj) = 
y^exp (—Tp(cj — 2co’o)^/2) and EPMF: = 

Apiujs uJi)Q{uJsi^i)‘ Within this definition the pulse 
duration Tp is related to value at full width half maxi¬ 
mum (FWHM) by Tp = 

In order to get a simplified expression for the bipho¬ 
ton wave function T(ci;s,ci;i) we approximate pump and 
SMF’s modes using Gaussian functions = 

Wpeyi-p(^-wlk\^, Up{ksT) = Wf exp , fi = 

respectively. The characteristic transverse mode di¬ 
ameters 2wp and 2wf are measured at 13.5% of in¬ 
tensity maximum. Next, the phase mismatch can be 
expanded up to first order in Tylor series in angular 
and spatial frequencies : Akz ~ Ps{^s ~ ^o) + 

Pi{ui - cJo) + dsx{ksx - ksxo) + dix{kix - kixo). For our 
particular setting the magnitudes of respective coeffi¬ 
cients are numerically equal within a good approxima¬ 
tion therefore we introduce the following simplifying as¬ 
sumption: dsx = —dix = dx and = “A = /^- This to¬ 


gether with the approximation for the sine function 10 
sinc(x) ~ exp(—x^/5) allows us to derive a simplified 
analytical formula for EPMF and consequently for the 
biphoton, which is again in a Gaussian form: 


T(ci;s, = A/'exp 


[ 24 
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with a characteristic width given by: 
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where JV stands for a scaling factor. Note that the wave 
function within this approximation is real. In the follow¬ 
ing part we juxtapose the approximate analytical model 
and precise numerical simulation results in order to dis¬ 
cuss conditions for optimal performance of a source. 

When a nonlinear crystal is pumped with a GW laser, 
the negative correlation of signal and idler photons’ fre¬ 
quencies originates in energy conservation relation. For 
pulsed pumping, the setup parameters can facilitate not 
only negative but also positive or no correlations. This, 
seemingly, energy conservation relation breaking, is a re¬ 
sult of a pump pulse spectral width and characteristics 
of down conversion process. 

In Fig. we present how positive and negative corre¬ 
lation can be achieved. Orientation of the biphoton wave 
function depends on phase matching condition and pump 
spectral amplitude. For photons with the same energy, 
EPMF is a function which maximum is approximately 
oriented in antidiagonal direction in frequency space as 



(c) L = 7.5 mm (d) L = 7.5 mm 
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Fig. 2. (a,c) The effective phase matching function 

0(co’s,co’i) and (b,d) the biphoton wavefunction T(ci;s,ci;i) 
for two exemplary crystal lengths (a,b) L = 1 mm and 
(c,d) L = 7.5 mm are depicted with contours. The func¬ 
tions are real. Green solid (dashed) lines represent pump 
pulse central frequency (spectral width). The simulations 
were performed for: = 50 fs, Wp = 100 /i m, 

Wf = 440 jam. 


depicted in Fig. |^a). A pump’s orientation is in the di¬ 
agonal direction, which is depicted by the green line and 
two green dashed lines representing pulse characteristic 
spectral width l/rp. The biphoton wave function is a 
product of those two as seen in Fig. [^b). 


We quantify the spectral correlation using the Pear¬ 
son coefficient, r, and taking the probability distribution 
based on the biphoton wave function Eq. 0: 


{{uJs — — <^o)} _ ^ '^P^+ /g-. 

- Wo)2)((wi - Wo)2) 1 + T^O-f 


We observe that a correlation is negative (positive) if the 
quantity a+Tp is > 1 (< 1). There is no correlation for 
cr+Tp = 1. The type depends on the relative widths of 
EPMF, cr+, and pulse spectral amplitude, l/rp . If the 
width of EPMF is larger (smaller) compared to pump’s 
as in Fig.j^a) (Fig.j^c)) the resulting biphoton spectral 
state is negatively (positively) correlated as in Fig. |^b) 
(Fig. [^d)). If both widths are equal, photons feature no 
correlation at all. 

Our analytical model shows that correlation param¬ 
eter r, introduced in Eq. ^ when combined with the 
definition Eq. 0, is a function of the pump pulse dura¬ 
tion Tp, the length of the crystal L and collected mode 


2 
























width Wf only. It does not depend on the pump spatial 
mode width Wp. Let us analyze the condition for spec¬ 
trally uncorrelated photons ^ ^)Tp = 1. It 

splits the space of setup parameters into regions of pos¬ 
itive and negative correlation. The formula shows that 
the thinner the crystal or the narrower collected mode 
profile, the shorter pulse is required for spectrally un¬ 
correlated photons. This observation also agrees with 
our numerical simulation results shown in Fig. We 
consider three exemplary settings involving two crystals 
thicknesses L = 7.5 mm and 2.5 mm and two pules dura¬ 
tions = 150 fs, 50 fs. The correlation parameter, 

r is depicted using contours, which are almost parallel 
to Wp axis for broad range of spatial mode profiles pa¬ 
rameters. 

It is favorable for some applications when the SMF 
coupled photon pair exhibits significant positive spec¬ 
tral correlation. In order to achieve that, EPMF width 
( 7 + has to be small. One way of doing that is by us¬ 
ing thick crystals as shown in Fig.[^b,d). The analytical 
approximation Eq. (§ shows a dependence on collected 
mode width 2wf and crystal thickness L. In the limit 
of long crystals or large SME collected mode profile, the 
EPME’s width cr+ is always finite. Therefore there is 
a broad range of experimentally available parameters, 
where cr+ is optimal for positive spectral correlation. It 
can be seen in Eig. |^a) that typical pulse duration of 
150 fs and thick crystal 7.5mm allows for positive corre¬ 
lation as high as r ^ 0.3 are possible. We observe that 
using longer crystals L > 7.5 mm with all other parame¬ 
ters fixed does not increase much the value of correlation 
parameter. On the other hand applying shorter pulse du¬ 
ration of 50 fs allows one to get r ^ 0.8 (0.2) for L = 7.5 
mm (2.5mm) , see Eig. ib, c). The use of shorter crys¬ 
tal is also beneficial in terms of SME coupling efficiency, 
which we will analyze later. However it limits the mag¬ 
nitude of spectral correlation. 

As a side note, we point out that the spectra of sig¬ 
nal and idler photons depend on a pules duration Tp as 
well as EPME’s width cr+. The widths, which are ap¬ 
proximately equal to each other, can be derived based 

on Eq. (1) as 




2dl 


+ 


/32L2 


^. It is easy to see 


that conoitions favorable for positive spectral correlation 
make the photons spectrally narrower for fixed pulse du¬ 
ration. 

Now we move on to an analysis of spectral correlation 
in conjunction with other source’s characteristics, which 
define its usefulnesses for certain applications. Those are: 
a SME coupled pair production rate Rc and a SME cou¬ 
pling efficiency 77 = Rc/VRsRi, where Rs and Ri are 
signal and idler photon production rates, respectively. 
We numerically analyze how those source char¬ 

acteristics depend on spatial mode widths of the pump 
Wp and collected mode Wf. Eig. shows our numerical 
simulation results: surface illustrates efficiency r] depen¬ 
dence, black contours on the top of the efficiency surface 


(a) L = 7.5 mm, Tp = 150 fs 
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Eig. 3. The SME coupling efficiency r] (surface plot), 
SME coupled pair production rate Rc (color scaling) and 
correlation parameter r (contours) as a function of pump 
2wp and collected 2wf spatial mode diameters. 


depict a correlation r and color scaling is for a SME pair 
production rate \og{Rc) . 

The SME coupled pair production rate Rc depends on 
spatial modes of the pump Wp and collected modes Wf. 
It increases very quickly with decreasing spatial mode 
widths, see red region in Eig. |^a,b,c) for small Wf and 
Wp. However it is almost flat for a broad range of pump 
and collected mode widths (yellow region). We also ob¬ 
serve it changes slowly with increasing crystal thickness 
L. This can be understood as a result of two competing 
effects: the longer the crystal, the higher the probability 
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to convert a pump photon into a pair, but on the other 
hand, the stronger the phase matching requirement lim¬ 
iting the spectral range of photons. 

Next, the SMF coupling efficiency, 77 , strongly depends 
on geometrical parameters, but exhibits weak depen¬ 
dence on spectral width of the pump l/r^. The results 
of our numerical simulations suggest that the highest 
coupling can be achieved when spatial modes of the 
pump and collected photons are both large and approxi¬ 
mately the same. For L = 7.5 mm (2.5 mm) and spatial 
modes diameters 2wp = 500/im 2wf = 500/im (340/im) 
it reaches r] ^ 0.6 (0.76) as can be seen in Fig. ||:a,c). We 
observe that the coupling efficiency, 77 , decreases when 
crystal thickness L gets bigger. This can be attributed 
to a transverse walk-off effect in birefringent type II crys¬ 
tal. 

One can use Fig. as a recipe for a source, which gen¬ 
erates pairs of a given correlation parameter r. Let us 
assume that the crystal thickness and pulse duration are 
fixed. The collected mode spatial profile is determined by 
a chosen correlation parameter. Next the optimal spatial 
pump mode can be chosen in order to achieve the best 
possible coupling efficiency or the SMF coupled pair pro¬ 
duction rate, Rc. Note also that Rc does change slowly 
in a wide range of spatial mode widths Wp and Wf^ which 
is in favor of the described optimization method. 

Summarizing, we discussed optimal settings for a 
source capable of generating photon pairs in telecom 
range featuring any kind of spectral correlation. For a 
given crystal thickness it is possible to control the spec¬ 
tral correlation by adjusting the pump pulse or collected 
modes characteristics. The correlation was analyzed to¬ 
gether with the SMF coupled pair production rate and 
coupling efficiency. Our numerical simulations showed 
that there is an asymptotic limit for the length of the 
crystal yielding maximal pair production rate. We also 
presented a recipe allowing one to simply choose opti¬ 
mal source parameters for a given correlation parame¬ 
ter and optimal both SMF coupled pair production rate 
and SMF coupling efficiency. In particular, the analyzed 
source can generate pairs featuring positive correlation. 
In addition to that the spectra of each individual photons 
got narrower with increasing positive correlation value 
for a fixed pulse duration. These two facts make it use¬ 
ful in applications where dispersion effects matter such 
as long distance fiber based quantum communication. 
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